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Abstract 

The spatial symmetry property of truncated birth-death processes studied in Di Crescenzo 
|S] is extended to a wider family of continuous-time Markov chains. We show that it yields 
simple expressions for first-passage-time densities and avoiding transition probabilities, and 
apply it to a bilateral birth-death process with jumps. It is finally proved that this symmetry 
property is preserved within the family of strongly similar Markov chains. 
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1 Introduction 

A spatial symmetry for the transition probabilities of truncated birth-death processes has been 
studied in Di Crescenzo jS] ■ Such a property leads to simple expressions for certain first-passage- 
time densities and avoiding transition probabilities. In this paper we aim to extend those results 
to a wider class of continuous-time Markov chains. 

Given a set {x n } of positive real numbers and the transition probabilities Pk,n(t) of a 
continuous-time Markov chain whose state-space is {0, 1, . . . , N} or Z, in Section 2 we introduce 
the following spatial symmetry property: 

PN-k,N-n{t) = —Pk,n(t)- (1) 
Xk 

In section 3 we point out some properties of first-passage-time densities and avoiding transition 
probabilities for Markov chains that are symmetric in the sense of (fl|). These properties allow 
one to obtain simple expressions for first-passage-time densities in terms of probability current 
functions, and for avoiding transition probabilities in terms of the 'free' transition probabilities. 
In Section 4 we then apply these results to a special bilateral birth-death process with jumps. 
Finally, in Section 5 we refer to the notion of strong similarity between the transition probabilities 
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of Markov chains, expressed by pk,n(t) = {fin I fik) Pk,n(t) (see Pollett [UUi and references therein) 
and show the following preservation result: if Pk,n{t) possesses the symmetry property then 
also Pk,n(t) does it. 



2 Symmetric Markov chains 

Let {X(t), t > 0} be a homogeneous continuous-time Markov chain on a state-space S. We shall 
assume that S = {0, 1, ... , N}, where iV is a fixed positive integer, or S = Z = {. . . , —1, 0, 1, . . .}. 
Let 

Pk,n{t) = Pr{X{r + t) =n\ X(t) = k}, k,n £ S; t,T>0 (2) 
be the stationary transition probabilities of X(t), satisfying the initial conditions 



Pk,n{0) = 5k,n 



1, k = n, 
0, k^n. 



Let Q be the infinitesimal generator of the transition function ©, i.e. the matrix whose (k, n)-th 
finite entries are: 

Qk,n = -^Pk,n{t) t=Q , (3) 

satisfying the following relations: (a) qk, n > for all k,n £ S such that k / n, (b) q n , n < for 
all n E S, and (c) J2 n eS ^ k ^ n = ® ^ or an ^ ^ 5. 

The spatial symmetry of Markov processes allows one to approach effectively the first- 
passage-time problem. Indeed, it has been often exploited by various authors to obtain closed- 
form results for first-passage-time distributions; see Giorno et al. |13j and Di Crescenzo et al. 
|H] for one-dimensional diffusion processes, Di Crescenzo et al. 7 for two-dimensional diffu- 
sion processes, and Di Crescenzo p] for a class of two-dimensional random walks. Moreover, in 
Di Crescenzo [H] a symmetry for truncated birth-death processes was expressed as in (jl]). with 
Xi suitably depending on the birth and death rates. Such symmetry notion can be extended to 
the wider class of continuous-time Markov chains considered above. Indeed, for a set of positive 
real numbers {x n ; n E 5} there holds: 

PN-k,N-n(t) = —Pk,n{t) for all k, n E S and t > (4) 

Xk 

if and only if 

qN-k,N-n = — Qk,n for all k,n£S. (5) 

Xk 

The proof is similar to that of Theorem 2.1 in Di Crescenzo and thus is omitted. 

Eq. (JIJ) focuses on a symmetry with respect to N/2, which identifies with the mid point of 
S when S = {0, 1, . . . , iV}. For each sample-path of X(t) from k to n there is a symmetric path 
from N — k to N — n, and the ratio of their probabilities is time-independent. Hence, in the 
following we shall say that X(t) possesses a central symmetry if relation @ is satisfied. 

Remark 2.1 If X(t) possesses a central symmetry, then 

Xn XN—k 



Xk X N . 



for all k,n E S. 
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An example of a Markov chain with finite state-space and a central symmetry is given hereafter. 



Example 2.1 Let X(t) be a continuous-time Markov chain with state-space S = {0,1,2,3}, 
with and 3 absorbing states, and infinitesimal generator 



Q 





ago + [3 -a(go + g 2 ) - (3 g 
a (3 






f3g ag 2 

-a(eo + g 2 ) - (3 go (ag + (3)g 




with a,(3,g > and go = 1 + Q. Then, X(t) has a central symmetry, with pjv- 
Q k ~ n Pk,n(t) for all k,n £ S and t > 0, and q^^k,N-n{t) = g k ~ n qk,n for all k,n £ S. 



■k,N—n 



Remark 2.2 If X(t) has a central symmetry and possesses a stationary distribution {7r ra , n £ 
S}, with lim pkn(t) = vr n > /or a// k,n € S, then the following statements hold: 

(a) Sequence {x n } is constant, so that pN_k,N-n{t) = Pk,n(t) f or a ^ k,n £ S and t > 0. 

(b) The stationary distribution is symmetric with respect to N/2, i.e. 

ftN-n = TTtj f or a tt n £ S. 

(c) Let X*(t) be the reversed process of X(t), obtained from X(t) when time is reversed, and 
characterized by rates and transition probabilities 



1k,n ~ Qn,ki 



7T r , 



Pn,k(t), 



k,n£S, t>0. 



Then, also X*(t) has a central symmetry, with p* N _ k N_ n (t) = P% n (t) for all k,n S S and t > 0. 
(d) Let D = {dk t n} be the deviation matrix of X{t), with elements (see Coolen-Schrijner and 
Van Doom 



dk,n 



\jPk,n{t) - 7T n ] dt, 



k,n E S. 



Then, D has a central symmetry, i.e. djv-fe Af-n = dk,n f or a ^ k,n£ S. 

An example of a Markov chain satisfying the assumptions of Remark 12.21 is the birth-death 
process on S with birth rate X n = a (N — n) and death rate fi n = an (see Giorno et al. |llj , or 
Section 4.1 of Di Crescenzo El). 



3 First-passage-time densities 

In this section we shall focus on the first-passage-time problem for Markov chains X(t) that 
have a central symmetry and that satisfy the following assumptions: 

(i) N = 2s, with s a positive integer; 

(ii) q id = q jti = 0, E ie 5_ %s > 0, T, je S+ 2?> > °> E iG 5_ Qs,i > and T, j£ S+ > for a11 
i G S- and j £ <S + , where 

S- = {n £ S; n < s}, S + = {n £ S; n > s}; 
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(in other words, if states i and j are separed by s then all sample-paths of X(t) from i to j, or 
from j to i, must cross s); 

(hi) the subchains defined on S- and S + are irreducibles. 

In addition, we introduce the following non-negative random variables: 



" l,S 



upward first-passage time of X{t) from state i £ S- to state s, 



T- s = downward first-passage time of X(t) from state j £ S+ to state s. 

We shall denote by gf s {t) and gj s (t) the corresponding probability density functions. Due to 
assumptions (i)-(iii), for all t > such densities satisfy the following renewal equations: 

(6) 
(7) 

(8) 
(9) 



Pi At) = / 9W)PsAt - tf) d#, i e 5_, je{ s }us + , 
io 

= / 9j,aW)P8,i(t - ■&) dtf, i E S—U {s}, j £ S , +. 



For all i > and k £ S let us now introduce the probability currents 

h+ s (t) = lim - P{X(t + t) = s, X(t) < s | X(0) = k}= V p fc>i (t) ?i>a> 

' tJ.0 T — ' 
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^ s (t) = lim - P{X(t + t) = s, X(t) > s | X(0) = k}= PkM Qj,s- 

They represent respectively the upward and downward entrance probability fluxes at state s at 
time t. Due to assumptions (i)-(iii) and Eqs. for i £ <S_, j £ 5+ and t > they satisfy 

the following integral equations: 

(10) 
(11) 



K>(t)= / gtmh- s (t-$)M, 

Jo 



Hereafter we extend Proposition 2.2 of Di Crescenzo |H] to the case of Markov chains. 

Proposition 3.1 Under assumptions (i)-(iii), for all i £ S— } j £ 5+ and t > the following 
equations hold: 



9ls(t) = hj t9 (t) - J q g- s ($) h- s (t - 0) dtf. 
Proof. For all t > and i £ 5_, making use of assumptions (i)-(iii) and Eq. (JSJ) we have 

^«.«(*) = E Pi,n{t) q n< s = hf s (t) + PiA 1 ) 9n,s- 



(12) 
(13) 



S 



nG{s}U5+ 



Hence, recalling (|SJ) we obtain 



dt 



91 



(*) + /%£(*) 

Jo 



E 

nG{s}US+ 



— p S)S (t - i9) - ^ p S)n (t - ??) q n . 

n£{s}US+ 



<&(0)p s , ft (i-0)d0 
dtf, 



(14) 
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where use of initial condition p ss (0) = 1 has been made. From Chapman-Kolmogorov forward 
equation we have 

n£{s}US+ 

so that Eq. (|14|) gives (|12j) . The proof of (|13|) goes along similar lines. ■ 

With reference to a Markov chain that has a central symmetry, we now come to the main 
result of this paper, expressing the first-passage-time densities through the symmetry state s as 
difference of probability currents © and @. 

Theorem 3.1 For a Markov chain that has a central symmetry and satisfies assumptions (i)- 
(iii), for all t > and k £ S there results: 

^, s W = 7iW' ( 15 ) 

Xk 

Moreover, for alii G 5_ ; j £ S + andt > the upward and downward first-passage-time densities 
through state s are given by 

= ht )S (t) - K s (t), gT s (t) = h~ s (t) - hf s (t). (16) 

Proof. Recalling that N = 2s, for t > we have 

h 2s-k,s( t )= P2s-k,j(t) qj, s (from©) 

= P2s-k,2s-i(t) q2s-i,s (setting j = 2s- i) 



ieS. 

X s 



%» ( from © and ©) 



■^k . c 

= ^Kjt). (from©) 

Xk 

Eq. (|15j) then holds. In particular, for k = s it implies that h~ s (t — i?) = hf s (t — 19) for all t > 
Hence, relations (|T6|) follow from Eqs. ([T0|) - (fT3|) , ■ 

For a Markov chain X(t) satisfying assumptions (i)-(iii) let us now introduce the s-avoiding 
transition probabilities: 

Pkl® = p W) = n > x (v) + s for a11 19 G (°> ') I x (°) = k } > 

where k,n £ S- U S + . We note that p^ n (t) is related to Pk, n (t) by 



(17) 



Pfe,n(i) - / 9 ks {^)Ps,n(t - #)dtf, k,n£S + . 



In the following theorem, for symmetric Markov chains two different expressions are given for 
pl s L(i) in terms of Pkn(t)- It extends Theorem 2.4 of Di Crescenzo |BJ; the proof is similar and 
therefore is omitted. 
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Theorem 3.2 Under the assumptions of Theorem V3.ll for t > and for k,n € 5_ U 5+ there 
holds: 

Pkliti) = Pk,n{ t ) ~ —P2s-k,n( t ) 
= Pk,n(t) ~ —Pk,2s-n(t)- 

We conclude this section by pointing out that for a Markov chain having a central symmetry, 
for all t > the following relations hold: 



4 A bilateral birth-death process with jumps 

In this section we shall apply the above results to a special symmetric Markov chain X(t) with 
state-space Z, characterized by the following transitions: (a) from n £ Z to n + 1 with rate 
A, (b) from n £ Z to n — 1 with rate /x, and (c) from n G Z — {0} to with rate a. Hence, 
X{t) is a bilateral birth-death process that includes jumps toward state 0. In order to obtain an 
expression for the transition probabilities Pk,n(t), we note that for all t > the following system 
holds: 

^PkA 1 ) = -(^ + A t + a )Pfc,n( t ) + Apfc in _i(i) +HPk,n+i(t), n G Z- {0}, 
^P*,o(*) = -(A + M)Pfc,o(*) +Apjfe_i(t) +fJ-Pk,i(t) +a^p iir (t), 

The probability generating function 

n=— oo 

is thus solution of 

0f 



— £T(z,t) = u( Z )#CM) +a, (19) 



where u(z) = —(A + [i + a) + \z + —, with initial condition 0) = z k . The unique solution 
of CU) is 

H(z,t) = H{z,0)e u{z)t + a [ e u ^ T dr. (20) 

Jo 

Hence, recalling that 

^- N /J nr. rv~i 



for 7 = 2y/\[i and /5 = y/X/fi, from 1)20(1 we obtain: 

+oo +oo 



JET(s,t) = e -( A +^+ Q )* ^ I n _ k (jt) (3 n ~ k s n + a 6 -(M-m+«)t ^ / n (7T)(/3s) n dr, (21) 



71= — OO 
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where I n (x) denotes the modified Bessel function of the first kind. Equating the coefficients of 
z n on both sides of (|21|) finally yields the transition probabilities 



n — k 



I n -k 2V A//t e 



-(\+H+a)t 



+ a 



f e -(A+M+a)T Jri (2V^ur)dr. (22 
Jo 



Note that (|2*2*)) can be expressed as 

Ph,n(t) = e" a *Pfc,n(i) + « / e~ aT p ,n(T)dT, 

Jo 

where, for all t > and k, n G Z, 



Pk,n(t) 



n — k 



/ n _ fc (2V^t) e -( A +^, 



(23) 



(24) 



is the transition probability of the Poisson bilateral birth-death process with birth rate A and 
death rate [i (see, for instance, Section 2.1 of Conolly [J). Assuming that the stationary prob- 
abilities 7r„ = lim pkn(t) exist for all n E Z, from ijTUl) we have 

t-^+oo ' 



E 



7T„ 



lim H(z,t) 



a 



n 



A(z - «i)(«2 - z) 



E 



where 



Hence, 



Zl,2 



A(^2 - «l) 

A + fj, + a ± V(A + + a) 2 - 4A^ 
2A 



Z \n 



E 

ra=0 



Z \n 
Z'2 



< z\ < 1 < z 2 - 



7T« 



A(^2 - Zl) 



az 2 



for n 



-1,-2 



1 1 ' ' ' 1 



(25) 



for n = 0, 1, 2, 



k A(z 2 - Zl) 

It is not hard to see that if A = \x then X(t) has a central symmetry with respect to state 0, 
with Xk = 1 for all k: 

P-k,-n{t) = Pk,n(t), Q-k-n = Qk,n 

for all t > and i,u6Z. Note that if A = /i, then z\ and 22 are reciprocal zeroes of u{z), so that 
the stationary distribution (f23|) is symmetric, i.e. ir n = 7r_ n for all n £ Z. Since = g^j = 0, 
g^o > and <7j,o > for all i,j € Z such that i < < j, and (?q,-i > and go,i > 0, this 
Markov chain satisfies assumptions (i)-(iii) for which is a symmetry state. In this case the 
first-passage-time densities through can be obtained via Theorem 13. II Indeed, if A = fi, making 
use of 112 2 H and of property I n (x) = I- n (x), for all t > and k = 1, 2, ... we have: 

+00 -1 
9k,o(t) = kfc,o(*) ~ h to(fi = E^J^) %'»o - E Pm(') «,o 

j'=1 i=— 00 
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Figure 1: On the left-hand are the plots of the downward first-passage-time density (I26j) for 
k = 3, A = 1 and a = 0.1,0.2,0.3, from bottom to top near the origin. On the right the 
0-avoiding transition probabilities 1)27(1 for k = 3, n = 1, A = 1 and a = 0.1,0.2,0.5, 1 (top to 
bottom) are indicated. 



-Pk,-i(t)] +a 



l= — oo 
+oo 



= e _ (2A+a)t | A [ Jjfe _ l( 2At) -J fc+1 (2At)] +aXj[/ Jfe _ i (2At) - I fc+i (2A t)] J. 
Furthermore, recalling (|18j) . in this special case for all t > and A; = 1,2,... there holds: 

In analogy with Theorem 13 . 21 and by virtue of (|22|). when A = /i, we have 

=Pk,n( t ) -P-kA 1 ) 



(26) 



Note that 



-(2A+a)* [j n _ fc ( 2 At) - 7 n+fc (2At)] , t > 0. 



(27) 

(28) 
(29) 



where is the transition probability of X(t) when A = fx. Functions (|26|) and (|27|) are 

shown in Figure 1 for some choices of the involved parameters. 

We finally remark that Eqs. (|23j) and (|29|) are in agreement with similar results for birth- 
death processes with catastrophes obtained in Di Crescenzo et al. [5] and |10j . 



5 Strong similarity 

The notion of similarity between stochastic processes has attracted the attention of several 
authors (see Giorno et al. J2]j f° r time-homogeneous diffusion processes, Gutierrez Jaimez et 
al. ^3] for time-nonhomogeneous diffusion processes, Di Crescenzo 0)0) and Lenin et al. [To] , 
for birth-death processes, and Pollett ^B], for Markov chains). Two continuous-time Markov 
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chains X(t) and X(t), with state-space S, are said to be strongly similar if their transition 
probabilities satisfy 

Pk,n{t) = Pjfc,n(*), for all t > and k,n eS, (30) 

Pk 

where {/3 n , n € S} is a suitable sequence of real positive numbers (we refer the reader to Pollett 
|16j . for further details). In the following theorem we state that if a Markov chain has a central 
symmetry, then any of its similar chains has a central symmetry as well. 

Theorem 5.1 LetX(t) andX{t) be strongly similar continuous-time Markov chains with state- 
space S; if X(t) has a central symmetry, then for allt>0 and k,n £ S one has: 

PN-k,N-n(t) = -^-Pk,n{t) 

Xk 

with 

x n — x n , n fc o. 

Pn 

The proof is an immediate consequence of assumed symmetry and similarity properties. 

Hereafter we show an application of Theorem 15.11 to a birth-death process having constant 
rates and state-space Z. 

Example 5.1 Let X(t) be the bilateral birth-death process with birth and death rates A and fi, 
respectively. From transition probabilities ()24j) it is not hard to see X{t) has a central symmetry 
with respect to 0, i.e. for all t > and k,n £ Z there results 

P-k-n(t) = —Pk,n(t), With X n = (- 

X k \fl 

The Markov chains that are strongly similar to X(t) constitute a family of bilateral birth- 
death processes characterized by birth and death rates (see Section 4 of Di Crescenzo and 
Example 3 of Pollett HHJ) 

A n = — 5— A, /i„ = — — ^, n e Z, 

Pn Pn 

and by transition probabilities (|3U|) . with Pk, n {~k) given in (|24j) and 

/ A\ n 

0n = 1 + V [-J , ri£Z, 

for all i] > 0. Due to Theorem l5.11 the family of strongly similar processes has a central symmetry 
with respect to 0: 

Xn~ u , ~ 0-n 1 + (J) " (\ 



p-k- n (t) = —Pk, n (t), with x n = — — x n = v \ n - , n G Z. 

Xk Pn 1 +77 ( AJ \fi; 
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